Upper limit on spontaneous supercurrents in Sr 2 Ru04 
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It is widely believed that the perovskite S^RuCU is an unconventional superconductor with broken 
time reversal symmetry. It has been predicted that superconductors with broken time reversal 
symmetry should have spontaneously generated supercurrents at edges and domain walls. We have 
done careful imaging of the magnetic fields above Sr2Ru04 single crystals using scanning Hall bar 
and SQUID microscopies, and see no evidence for such spontaneously generated supercurrents. We 
use the results from our magnetic imaging to place upper limits on the spontaneously generated 
supercurrents at edges and domain walls as a function of domain size. For a single domain, this 
upper limit is below the predicted signal by two orders of magnitude. We speculate on the causes and 
implications of the lack of large spontaneous supercurrents in this very interesting superconducting 
system. 



I. INTRODUCTION 



The perovskite superconductor Sr 2 Ru0^i (T C =1.5K) 
is believed for a number of reasons to have uncon- 
ventional pairing symmetry^ Muon spin resonance ex- 
periments are consistent with the generation of large 
but sparse internal magnetic fields when Sr2Ru04 be- 
comes superconducting, indicating a superconducting 
state with broken time reversal symmetry.— This finding 
is supported by the observation of the onset at the super- 
conducting transition temperature of a Kerr effect rota- 
tion of light polarization upon reflection^ consistent with 
large (~50-100/im) domains with broken chiral symme- 
try. The first phase sensitive Josephson tunneling mea- 
surements suggested a static 7r-phase shift between oppo- 
site faces of a Sr 2 Ru04 single crystal^ a result obtained 
only when the sample was prepared by controlled slow 
cooling, implying that the sample could be made to pos- 
sess a small even number of domain walls separating the 
opposite faces. On the other hand, more recent Joseph- 
sen tunneling measurements point toward small (~l/im) 
dynamic order parameter domains^ Magnetic imaging of 
the ab face of a single crystal of Sr 2 Ru04 using a micron 
sized SQUID shows vortex coalescence on a scale of ^10 
/im that may be related to a domain structure in the su- 
perconducting order parameter! 7 - These observations, as 
well as others, have been interpreted in terms of a su- 
perconducting order parameter with spin triplet, 8 chiral 
p x ± ip y Cooper pairing symmetry. However, Hall bar 
microscopy measurements^ did not observe the magnetic 
fields expected at the surface or edges of a superconduc- 
tor with broken chiral symmetry. Here, we report on 
new scanning SQUID microscope measurements, and we 
also further analyze some of the Hall bar measurements, 



showing that, if they exist, the local fields at the surfaces 
and edges of Sr 2 Ru04 single crystals are much smaller 
than those expected at the surfaces and edges of a chiral 
p-wave superconductor. 

Although p x ± ip y pairing symmetry is fully gapped, 
specific heat,— nuclear relaxation rate,— and thermal 
conductivity measurement a 12 i 13 ' all show a power law 
temperature dependence, suggesting the presence of line 
nodes. Among other suggestions, one possibility is that 
the 7 band, the band with the primary contribution to 
superconductivity in Sr 2 Ru04, has nodeless p x ±ip y pair- 
ing symmetry, but induces superconductivity with a line 
of nodes in the other (a and j3) bandsJ^ Anisotropy in 
the gap function— has been supported experimentally by 
specific hea t 17 i 18 and ultrasound attenuation^ measure- 
ments. 

The issue of the broken time reversal symmetry in the 
superconducting state of Sr 2 Ru04, aside from intrinsic 
interest, has taken on new urgency with several proposals 
for error tolerant quantum logic elements taking advan- 
tage of this propertyi 20 ' 21 i 22 

There are useful analogies between a chiral p x ± ip y 
superconductor and a ferromagnetics 3 - A single domain 
ferromagnet has a uniform magnetization which is equiv- 
alent to the field produced by a current sheet circulating 
around the surface, in the appropriate geometry, while 
a single domain p x + ip y superconductor carries an in- 
trinsic angular momentum of fi per Cooper pairjSi which 
one would expect to lead to an actual surface current 
sheet, confined within a healing length proportional to 
the coherence length of the surfaced However, the field 
generated by this current must be screened inside the 
superconductor by a diamagnetic shielding current flow- 
ing within the penetration depth of the surface, so that 
B = inside the superconductor. The net result is 
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a spontaneous magnetization within the healing length 
plus the penetration depth of the sample edges, which 
is greatly reduced from that expected from the simple 
ferromagnetic analogy, but which is still substantial and, 
using parameters appropriate to Sr2Ru04, predicted to 
give rise to local fields as large as lmT under certain 
assumptions! 26 ! 27 The superconductor can also support 
domains in which regions of p x + ip y coexist with regions 
of p x — ip y order. Although the net magnetization van- 
ishes at the boundaries between such domains, the local 
fields, which extend over the penetration depth on ci- 
ther side of the wall, can be as large as 2mTi 26 i 27 Other 
than direct phase sensitive measurements, the detection 
of such fields would be one of most direct confirmations 
of a superconducting order parameter with time reversal 
symmetry breaking since the spontaneous boundary and 
domain wall supercurrents are expected by symmetry 2 ^ 
and would have no other obvious explanation. It is there- 
fore appropriate to attempt to image the magnetic fields 
arising from these spontaneous supercurrents using scan- 
ning magnetic microscopy. 



II. MAGNETIC IMAGING 

We have performed scanning magnetic imaging of the 
ab and ac faces of single crystals of Sr2Ru04. The mag- 
netic images reported here were made at Stanford with a 
dilution refrigerator based Hall bar/SQUID microscope 29 
with a base temperature below lOOmK, and at IBM 
with a 3 He based scanning SQUID microscope with a 
base temperature below 300mK. Our SQUID sensors had 
square pickup loops 8 /im on a side; the Hall bars had 
roughly square effective areas 0.5 /im on a side. The Hall 
bar measurements were made in a residual field of about 
2.5/iT; the SQUID measurements were made in a resid- 
ual field of 75nT, compensated for fields perpendicular 
to the scanning direction to less than lOnT using a small 
Helmholtz coil. The SQUID measurements were made 
after cooling the samples through the superconducting 
transition temperature at a rate of about lmK/sec. Some 
of the Hall bar data discussed in this paper has been re- 
ported previously*?- However, here we make a more quan- 
titative comparison of this data with theory. 

The Si'2Ru04 single crystals used in our experiments 
were grown using a floating zone method^ The samples 
used for the IBM SQUID measurements were mounted 
in epoxy and polished so that either the ab or ac face 
was part of a smooth plane, allowing scanning across 
the edges of the crystal^ Some of the samples used for 
SQUID microscopy were the same as for phase sensi- 
tive experiments on the pairing symmetry of Sr2Ru04, 5 
and had layers of SiO and Auo.slno.s (T c = 0.4-0.5K) de- 
posited on some of the crystal faces perpendicular to the 
scanned face. These additional layers should have had 
no effect on the magnetic imaging experiments reported 
here. The critical temperature of the crystals were mea- 
sured to be > 1.4K using scanning and bulk susceptom- 
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FIG. 1: SQUID microscope image of the ab face of a S^RuCU 
single crystal, cooled in a field B z <10 nT and imaged at 
T=0.27K with an 8 (jm square pickup loop, a Pseudocolor 
image with full-scale variation of 0.2 "l>o ($o = h/2e) in mag- 
netic flux through the SQUID pickup loop. The dashed line 
in a shows the outlines of the crystal, b Same image as a, but 
with the pseudocolor scale expanded to 0.02 $q. The dashed 
line in b shows the line traced by the cross-section in c. The 
dashed rectangle in b shows the area of the image expanded 
in d. e is a histogram of pixel values for the data displayed 
in d. 



etry measurements. 

Figure [T] shows a SQUID microscope image of the ab 
face of a Sr 2 Ru04 single crystal. The largest feature ev- 
ident in this image (Fig. [Jti) is an isolated Abrikosov 
vortex. When the pseudocolor scale is expanded to 
A$ s = 0.02<f>o (Fig. [lb) magnetic features become ap- 
parent in the epoxy and along the edges of the crystal. 
We believe these features are not due to the supercon- 
ductivity of the Sr 2 Ru04 because they are unchanged 
from cooldown to cooldown in different fields. Figure [2] 
compares images from 3 different cooldowns of the same 
crystal, in nominal ambient plus compensating fields per- 
pendicular to the scanning plane of zero (Fig. [2Kh 10nT 
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(Fig. EJd), and 15nT (Fig. [2J:). The number and po- 
sitions of the Abrikosov vortices in the top (ab) face 
of the crystal, and an interlayer vortex emerging from 
the left (ac face) edge (Fig. O)) of the crystal change 
from cooldown to cooldown, but the sharp features at the 
edge of the sample are remarkably reproducible. These 
edge features may be the result of the polishing process, 
such as topographical or magnetic features from particles 
trapped in the epoxy. Note that features very similar to 
the edge features are apparent in the epoxy far from the 
sample edge. Above the sample itself the flux image is 
relatively smooth, with a broad background (Fig. [1};). 
We believe that this broad background is the result of 
magnetic flux coupled into the SQUID through sections 
outside of the pickup loop. A clear demonstration of this 
effect appears in Ref. [3lJ. On top of the broad back- 
ground, two steps in the cross-section (Fig. [TJ:) corre- 
spond to the edges of the crystal. We believe that these 
steps are due to small supercurrents circulating around 
the entire sample due to uncompensated residual fields 
(see Figure [5]). Figure [T]i shows a magnified image of a 
section of the crystal (indicated by the box in Fig. [TJd), 
with no magnetic features larger than a few m<E>o over an 
area of several hundred microns on a side. 

Similar results were obtained when SQUID microscope 
images were taken of the ac face of a Sr2Ru04 single 
crystal (Figure [3]) . In this case there were a number of 
interlayer vortices with flux both emerging from and en- 
tering into the crystal surface near the left edge of the 
crystal (Fig. [3Jd). Just as for the ab face, there were sharp 
magnetic features along the edges of the crystal and in 
the epoxy which did not appear to be correlated with 
the superconductivity of the Sr2Ru04, as well as broad 
magnetic backgrounds, but sharp magnetic features were 
absent from large areas of the crystal face. 

The samples used in the Hall bar measurements were 
cleaved. 1/im diameter, ~1 /im deep holes were milled on 
a 20 /im grid on the upper surface using a focussed ion 
beam, to create artificial edges.— Figure H] shows a scan- 
ning Hall bar image of the ab face, with a regular array 
of 1 micron holes at a pitch of 20 microns, of a S^RuCU 
single crystal. There are a few Abrikosov vortices appar- 
ent in this image, but the area away from these vortices 
is featureless. In particular, no features were observed 
in connection with the edges or interiors of the 1 micron 
holes. Since the 1 /xm deep holes did not even act as effec- 
tive pinning centers for the vortices, they may not have 
served as significant singularities to create edge currents. 
The outer edges of the crystal were not scanned in the 
Hall bar measurements. 



III. MODELING 

Matsumoto and Sigrist 26 (MS) have solved the 
Bogoliubov-De Gennes equations using a quasi-classical 
approximation for the cases of an edge between a semi- 
infinite, ideal p x 4- ip y superconductor and vacuum, and a 
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FIG. 2: Comparison of SQUID microscope images of an ab 
face of a Sr2Ru04 crystal after three different cooldowns in 
slightly different magnetic fields. 



domain boundary between a,p x + ip y superconductor and 
a Px — ipy superconductor. Their solutions are fully self- 
consistent so that they include the effect of screening cur- 
rents. They predict substantial supercurrents and conse- 
quent magnetic fields spontaneously generated at edges 
and domain boundaries. For example, the peak mag- 
netic fields in these calculations correspond to 1 mT for 
edges, and 2 mT for domain walls using values for the co- 
herence length (£o=66 nm) and penetration depth (Xl — 
190 nm) suitable for Sr2Ru04. However, some modeling 



4 



500/mi 




-0.02 — 1 1 

-1 1 
x10 3 

Position (/jm) 



FIG. 3: SQUID microscope image of the ac face of a S^RuCU 
single crystal cooled in nominally zero field and imaged at 
T=0.27K with an 8 fiva square pickup loop, a Pseudocolor 
image A<E> S = 0.8$o- The dashed line in a shows the positions 
of the outer edges of the crystal, b Same image as a, but with 
A$ s = 0.08"3>o A few interlayer vortices with both positive 
and negative signs are visible near the lower-left edge of the 
crystal. The dashed line in b is along the a-axis and shows the 
data traced by the cross-section in c. The arrows in c indicate 
the edges of the crystal. The dashed square in b shows the 
area of the image expanded in d. The diagonal stripes visible 
in d are due to 60 Hz noise, e is a histogram of pixel values 
for the data displayed in d. 



is required to compare our experimental results with the 
MS predictions because we measure the magnetic fields 
above the surface, rather than inside the sample. 

The simplest approach to this problem is to assume 
that the magnetic fields at the surface of the sample are 
the same as those in the bulk. This neglects field spread- 
ing and any change in superconducting shielding due to 
the finite sample geometry. However, in our case the size 
of the magnetic sensor and its spacing from the sample 
are large relative to the coherence length and penetra- 
tion depth, so that the field averaging from these effects 
are larger than the additional effects of field spreading 
and changes in superconducting shielding. The field av- 
eraging effects from finite sensor size and height can be 
shown rigorously to be larger than field spreading, for ex- 
ample, in the similar problem of vortex fields spreading 
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FIG. 4: a Scanning Hall bar image of ab face of Sr2Ru04 
single crystal, cooled in ~ 2.5//T and imaged at a temperature 
below lOOmK using a Hall bar with a sensor area 0.5/xm on a 
side. In this image the mean of each scan line was subtracted 
from the raw data to remove slow drift in the sensor Hall 
voltage, b Same area as a but with an expanded pseudocolor 
scale. The dashed line in b shows the line traced by the data 
cross-section in c. The dashed square in b shows the area for 
which a histogram of pixel inductance values are displayed in 
d. 

from the surface of a superconductor and imaged with a 
SQUID microscopei 32 ! 33 In the remainder of this section, 
we will neglect changes in the currents near the surface 
due to the finite sample geometry. We will show below 
that the effect of finite sample geometry only leads to 
suppression of the expected signal by 30% compared to 
what is expected from the edge currents of an infinite 
sample. The finite sample geometry effects for edge cur- 
rents are expected to be similar to those for domain walls 
and are also discussed in the following section on surface 
screening effects. 

It is well known^i that if the normal component of 
the magnetic field B z (x, y, z) is known at all points of a 
surface z — the magnetic field in free space at a height 
z above that surface is given by 

B z (k Xl k yi z) = B z (k x ,k y ,z = 0)er kz , (1) 

where B z (k y ,k y ,z) is the 2-dimensional Fourier trans- 
form of B z (x,y,z) and k = ^Jk% + k y . To model the 

magnetic signals in our experimental SQUID and Hall 
bar microscope geometries, we assume a particular do- 
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FIG. 5: Cross-section through the image of the ab face of 
Sr2Cu04 displayed in Figure [T] (solid line). The short-dashed 
line is the prediction for a superconducting disk in a uniform 
residual field of 3 nT. The long-dashed line (with a peak at 
$ 8 /$o = 1.1) is the prediction for a single domain p x + ip y 
superconductor of the extended-Matsumoto-Sigrist model as 
described in the text, assuming a square pickup loop 8(j,m 
on a side, at a height of 3/xm above the sample. Here the 
superconductor is positioned to the left of with epoxy 

to the right. 
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main structure with the magnetic fields, B z , at each edge 
and domain boundary, at the surface z = 0, taken to be 
those predicted by Matsumoto and Sigrist^ for an infi- 
nite sample. We then propagate the fields to a height 
z using Eq. [TJ integrate over an area appropriate for 
the SQUID or Hall bar sensor to obtain a magnetic flux, 
and divide by the area of the sensor for the case of the 
Hall bar to get an average magnetic field. We will re- 
fer to this model as the "extended-Matsumoto-Sigrist" 
model to distinquish it both from the prediction made 
by Matsumoto-Sigrist for an ideal (infinite) geometry 
and from the more accurate model which includes ad- 
ditional screening effects due to the finite geometry, as 
discussed in the next section. The original Matsumoto- 
Sigrist results are scaled in field by B c — $ /2\/27r£oA.L, 
where $ = h/2e is the superconducting flux quantum, 
£o is the coherence length and Al is the London pene- 
tration depth. For the modeling presented here we take 
£o = 66nm and Xl — 190nm.— 

Figure [5] compares the results of this calculation (long- 
dashed line) with the experimental cross-section of the 
image shown in Figure [1] (solid line). Also shown for 
comparision is the predicted cross-section for an ideal 
superconducting disk in a uniform residual field 35 of 3 
nT. The small steps in flux at the edges of the crystals 
in Figures [T] and [3] can be attributed to shielding of a 
very small residual background field. These steps are 
much smaller than the peaks predicted by the extended- 
Matsumoto-Sigrist model for a single domain. 



FIG. 6: Predicted magnetic fluxes through an 8 /im square 
pickup loop, 3 fim above the sample surface, for a 64 /im 
square p x ± p y superconductor with various domain sizes, us- 
ing the predictions for the edge and domain wall currents of 
Matsumoto and Sigrist as described in the text. The dashed 
lines in the figure show the positions of the cross-sections dis- 
played in Fig. [7] 



Figure [6] shows the results from the modeling outlined 
above for a series of domain sizes using parameters ap- 
propriate for our SQUID measurements. In these cal- 
culations it was assumed that the domains were square, 
and extended infinitely far in the negative z direction 
(perpendicular to the crystal face). Fig. shows cross- 
sections through the modeling results as indicated by the 
dashed lines in Fig. O As expected, the magnetic fields 
above the edges and domain boundaries are averaged over 
a length set by both the height of the sensor above the 
sample surface and its size. This leads to a rapid decrease 
in the predicted signal when the domains become smaller 
than a critical length. (In this modeling the magnetic sig- 
nal for a domain size of 4 microns vanishes everywhere 
except at the sample corners because, due to the sym- 
metry of the domains with respect to the sensor, there 
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FIG. 7: Cross-sections through the modeling images of Fig 
[5] for various domain sizes. 



are exactly as many positive as negative contributions to 
the flux through the 8 /im diameter pickup loop. For this 
reason, we show the predicted flux for 32/6=5.3 rather 
than 4 micron domains.) The calculated peak values for 
the SQUID flux signal for edges and domain boundaries 
are plotted in Fig. [8^ as a function of domain size. The 
lower dashed line in Fig. [5£i is an estimate of the noise 
in the SQUID images above the interior of the crystals, 
taken to be the rms noise of the flux distribution shown 
in Fig. [3^ (2.5m$o)- The upper dashed line is the rms 
value of the flux distribution above the sample edges in 
Fig. [Th. (8.5m$ )- Comparable modeling results using 
parameters appropriate for our Hall bar measurements 
are shown in Fig. [5Jd. In this case the dashed line repre- 
sents the rms noise value of the field distribution in Fig. 
HJi (3.5/xT). We do not display an experimental limit on 
the possible edge currents set by the Hall bar experi- 
ments because of uncertainties associated with the hole 
geometry and surface damage induced by the focussed 
ion beam in these experiments^ 

In order to place limits on the possible field magnitude 
and domain sizes consistent with our results, we assume 
that the magnitude of the spontaneous supercurrents can 
vary, but that the spatial distribution of spontaneous su- 
percurrents is as calculated by Matsumoto and Sigrist. 
With this assumption we can scale the results, for exam- 
ple, in Figs. [5k,b vertically by assuming a scaling field 
B s different from B c = 4>n/2\/27r£oAL. In order for the 
spontaneous supercurrents to be unobservable in our ex- 
periments, the scaling factor and domain size must be 
in the region below and to the left of the lines in Figs. 
(5J;,d. Either the spontaneous currents are substantially 
smaller than calculated from the extended-Matsumoto- 
Sigrist model, or the domains are small. For example, 
for the SQUID measurements, the magnitude of the su- 
percurrents at the edge must be a factor of 100 smaller 
than those prediced by MS if the domains are 10 or more 
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FIG. 8: a Plots of the predicted peak flux signals for an 8/im 
square SQUID pickup loop, 3 fim above the sample surface, 
for a 64 /im square p x ± ip y superconductor with various do- 
main sizes, using the predictions of Matsumoto and Sigrist 26 
for the spontaneously generated edge and domain supercur- 
rents. The dashed lines represent the estimated SQUID noise 
in the measurements within the sample (lower line) and at 
the sample edges (upper line), b Plots of the predicted peak 
fields for a square Hall bar 0.5 (im on a side, 1.2pm above the 
sample surface, with the corresponding Hall bar noise floor, 
c Upper limits on the size of the scaling fields B s , normalized 
by B c = &o/2\/^Tr£ t o^L, as a function of domain size, given 
by our failure to observe spontaneously generated supercur- 
rents at edges and domain walls in the SQUID measurements. 
In this figure the extended-Matsumoto-Sigrist predictions are 
represented by B a /B c —1. d Upper limits on B s /B c as a func- 
tion of domain size set by the Hall bar measurements. 
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IV. SURFACE SCREENING EFFECTS 

In our modeling we have neglected the fact that the 
magnetic fields at the sample surface will be somewhat 
reduced from their bulk values. In principle, one can cal- 
culate the surface fields by self-consistently solving the 
Bogliubov-de Gennes equations in the appropriate ge- 
ometry. Here, we simply estimate the errors involved in 
neglecting surface screening effects using a London ap- 
proach. Following Ref. [3y, the superconductor is as- 
sumed to fill the half-space z < 0. If the change in the 
penetration depth close to the surface is neglected, the 
magnetic field B inside the superconductor can be de- 
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composed as B = Bo + Bi, where Bo is the particular 
solution given by Matsumoto and Sigrist^S of the inho- 
mogeneous London's equation for a domain wall and Bi 
is a general homogenous solution chosen to satisfy the 
matching conditions at z = 0. The London's equation 
for the particular solution can be written as^& 

k(K+k)<f> K = KeVB (k, 0)+ik- (B (k, 0) - 47rM(k, 

(2) 

where k = ^jk\ + fc 2 , K = \Jk 2 + 1/A 2 , the magnetic 

field B above the superconductor is given by B = —V$>k> 
Bo and M are the 2-dimensional Fourier transforms 
in x and y of the inhomogeneous solution to London's 
equation and the volume magnetization respectively, and 
k = k x e x + k y e y . However, if the domain walls are as- 
sumed to be parallel to the z axis, both Bo and M have 
only z components, and Eq. ([2]) reduces to 



K 



k(k + K) 



S O2 (k,0), 



where B 0z is the z-component of B . Then 



B(k,z) 



(ik + ke z 



k(k + K) 



K B 0z (k,0)e- kz . 



(3) 



(4) 



In our case we are only interested in the z-component of 
the field outside of the superconductor, which takes the 
particularly simple form 



B z (k,z) = 



K 



k + K 



B 0z (K0)e 



-kz 



(5) 



The modeling in the previous section, which neglects 
surface shielding, is equivalent to Eq. [5] in the limit 
A — > 0. Fig. [5] shows the effects of surface screening 
on the fields predicted for a single domain boundary for 
parameters appropriate for our Hall bar measurements. 
Even in this case the effects of screening are relatively 
small, because the penetration depth is smaller than the 
measuring height and size of the Hall bar. Surface screen- 
ing effects would be even smaller (a few percent) for the 
case of SQUID imaging, because of the larger size of the 
sensor. The geometry for considering the effects of super- 
conducting shielding on the edge fields is more complex 
than for the case of the domain boundary, as one needs 
to consider a superconductor bounded by both z and at 
least one of x or y. However, again the edge and surface 
effects will be confined on the scale of the penetration 
depth which is much smaller than the distance to the 
probe or the probe size. Therefore, we do not believe 
that the simple model presented above will be more than 
a factor of two different from a full calculation. 



V. DISCUSSION 

If the superconductivity of Sr2Ru04 breaks time- 
reversal symmetry, it should spontaneously generate su- 
percurrents at domain boundaries and sample edges. The 
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FIG. 9: Comparison of the predicted magnetic field sensed 
by a 0.5 /jm square Hall bar, 1.2/xm above a single domain 
in a p x ± ip y superconductor, using the spontaneous do- 
main currents predicted by Matsumoto and Sigrist,— with 
(dashed line) and without (solid line) surface screening ef- 
fects as described in the text. The solid line corresponds to 
the extended-Matsumoto-Sigrist model. 



fact that no magnetic fields due to such supercurrents 
were observed using scanning magnetic microscopy places 
significant limits on the size of these currents and the size 
of the domains, as shown in Fig. [5] In particular, from 
the combined Hall bar and SQUID measurements, we 
conclude that if the spontaneous supercurrents at a do- 
main wall are of the size expected from the calculations 
of Matsumoto and Sigrist 2 - and the modelling done here, 
one can set a conservative upper limit on the domain size 
of 1.5 microns for both interior and edge domains. Al- 
ternatively, if the domains intersecting the ab face are 
10 microns or more in size, we conclude that the sponta- 
neous supercurrents at edges are a factor of 100 smaller 
than expected from the calculations of Matsumoto and 
Sigrist combined with our modelling. 

Calculations of the self-consistent screening currents 
employed in our modelling have assumed an ideal p x ±ip y 
superconducting gap symmetry^ However, for S^RuOi, 
the gap in the ab plane is believed to be anisotropic ! 17 ! 18 
In addition, three different bands contribute to the Fermi 
surface in S^RuO^ These properties are likely to im- 
pact on the magnitude of the self-consistent screening 
currents, although, a priori, it is not clear whether the 
magnitude would be increased or decreased from the val- 
ues calculated by MS. On the other hand, muon spin 
resonance observed internal fields which are roughly con- 
sistent with the predicted values^ If these observed fields 
are due to internal domain walls, it suggests that the 
surface currents must be reduced by two or more orders 
of magnitude from their bulk values if the domains are 
larger than 10 microns. It is difficult to imagine what 
could so strongly reduce the surface fields at the ab sur- 
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face due to domain boundaries. The surface screening 
effects are small; the surfaces are cleaved; and roughness 
even to the depth of a hundred angstroms or so will not 
substantially reduce the fields detected at the Hall probe 
or SQUID. This suggests that domains intersecting the 
ab plane are either so sparse as to not have been scanned 
or are smaller than a few microns. Another possibility 
is that the domains intersecting the ab surface are shal- 
low, with a depth along the c-axis noticeably less than 
the penetration depth. In this case, the spontaneous cur- 
rents and fields could be too weak and spread out in the 
ab layers to be detected. However, we note that either 
small domains or domains shallower than the optical skin 
depth would also interfere with observations of the Kerr 
effect rotation. 

Sufficient roughness of the ac or be faces can be ex- 
pected to have a more noticeable effect on the edge or 
boundary currents. The samples used for the SQUID 
measurements on these faces were polished and AFM 
imaging on typical samples show them to be smooth to 
5 nm (rms)j^ MS assumed specular scattering from the 
edge in which case one component of the order param- 
eter is suppressed while the other component is slightly 
enhanced. For diffuse scattering from a rough edge, both 
components will be suppressed and this will reduce the 
surface currents and resulting magnetic fields. Although 
self-consistent calculations have not been carried out for 
this case, the effect of surface roughness on the two com- 
ponent order parameter has been studied^ and one finds 
that the two components heal over quite different length 
scales. Using Ginzburg-Landau and London theory to es- 
timate the resulting change in the surface magnetization, 
one finds that, even for completely diffuse scattering from 
a rough surface, the reduction in surface magnetization 
is less than 30% M 

Domain walls cost energy because they disrupt the su- 
perconducting order. Unlike a ferromagnet, there is no 
balancing of this energy due to dipolar forces because 
spontaneous screening currents ensure that the magnetic 
field, or local magnetization, is zero inside the supercon- 
ductor. Therefore, in principle, a single domain p x + ip y 
superconductor is possible. However, domains will nat- 
urally form as the sample is cooled through T c and as 
extended objects, these domains are susceptible to pin- 
ning by defects and impurities in the sample. Therefore, 
one expects domains to be present although their den- 
sity may be controlled by sample purity and slow cooling 
in a field. Muon spin resonance experiments were inter- 
preted as evidence for dilute domains^ Kerr effect mea- 
surements suggested domain sizes in the range of 50 to 
100 microns, - while the first phase sensitive Josephson 
tunneling measurements are consistent with no domains 
(or a small even number of domain walls between oppos- 
ing faces of the crystal) £ On the other hand, more recent 
Josephson tunneling measurements were interpreted as 
evidence for dynamic domains of ~1 micron on average,— 
although one would extract larger domain sizes if finite 
domains perpendicular to the c-axis were included in the 



modelling. All of these measurements, except for muon 
spin resonance, would see reduced signals if the domain 
size along the c-axis becomes small and this would af- 
fect the measurements reported here as well. Unless the 
fields at domain walls are reduced by more than an or- 
der of magnitude in size from the predicted values, the 
Hall bar measurements suggest domain sizes of either less 
than 1.5 microns in size over the ab face or large enough 
that no domain wall fell in the lOOx 100 /xm scan area. 

Earlier work has reported that large domains can be 
flipped by fields of the order of a mT or larger— and that 
small surface domains are influenced by fields < 0.1/iT. 6 
While the data presented here was taken on samples 
cooled in fields less than 2.5/xT, Hall data taken on sam- 
ples cooled in up to a mT was very similar to that shown 
here except for the presence of more trapped vortices.— 
In principle, very fast domain wall motion could result in 
zero time-averaged edge current and zero time-averaged 
domain wall current. However, previous experiments^^ 
suggest that the domain wall motion would be slow on 
our experimental time scale, which is 10 secconds per 
line scan for the scanning SQUID microscope data shown 
here. Therefore, it is unlikely that dynamic behavior of 
the domains prevents the observation of the signal in this 
experiment. 

In conclusion, scanning magnetic microscopy measure- 
ments place quite severe limits on the size of edge cur- 
rents and/or on domain sizes in S^RuO-j. The differ- 
ent experimental results taken as evidence for p x + ip y 
pairing come to quite different conclusions about domain 
sizes. Since there are now detailed predictions for the 
field profile in the vicinity of domain walls in the bulk, 
muon spin resonance could now, in principle, provide de- 
tailed information about the validity of these predictions 
as well as quantitative information about the density of 
domains in the bulk. In addition, either slow muons^ or 
beta-NMR^i could be used to probe the surface region 
and to look for fields due to spontaneous edge currents 
as well as domains near the surface. Scanning magnetic 
microscopy is still one of the most direct probes of do- 
mains intersecting the surface and of edge currents and 
further improvements in sensitivity may either confirm 
or rule out their existence. 
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